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Using 535 x 10 6 B -meson pairs collected by the Belle detector at the KEKB e + e collider, we 
measure branching fractions of (7.16 ± O.lO(stat) ± 0.60(syst)) x 10 -4 for B + —> J/4>K + tv + tv~ and 
(4.31 ± 0.20(stat) ± 0.50(syst)) x 10" 4 for B+ -s> ip' K + -k + -k~ . We perform amplitude analyses to 
determine the resonant structure of the K + -k + tt~ final state in B + — > J/i/jK + tt + -k~ and B + — ¥ 
ip'K + 7r + TT~ and find that the ifi(1270) is a prominent component of both decay modes. There 
is significant interference among the different intermediate states, which leads, in particular, to 
a striking distortion of the p line shape due to the uj. Based on the results of the fit to the 
B + -)■ J/tpK + TT + TT~ data, the relative decay fractions of the 7fi(1270) to Kp, Ku, and A"*(892)tt 
are consistent with previous measurements, but the decay fraction to ifo(1430) is significantly 
smaller. Finally, by floating the mass and width of the 7^1(1270) in an additional fit of the B + — > 
J/ipK + ir + Tr~ data, we measure a mass of (1248.1 ± 3.3(stat) ± 1.4(syst)) MeV/c 2 and a width of 
(119.5 ± 5.2(stat) ± 6.7(syst)) MeV/c 2 for the ifi(1270). 



PACS numbers: 13.25.Hw, 13.25.Es, 14.40.Df 

I. INTRODUCTION 

The large number of P-meson decays observed at B 
factories allows detailed studies of the intermediate-state 
resonances involved in these decays. This paper analyzes 
the structure of the K + ir + ir~ final state in the decays 
B+ -> J/ipK+n+n- and B+ -> tp' K+tt+tt-R Kaon ex- 
citations that decay to a Kirir final state are difficult 
to distinguish based on the mass of the Kitit system 
alone, owing to their overlapping line shapes0 In this 
analysis, data are therefore fitted in the three dimen- 
sions M 2 (Ktttv), M 2 (Ktt), and M 2 (irir), which are the 
squared invariant masses of the A + 7r + 7r~, K + tt~ and 
7r + 7r~ systems, respectively. An unbinned maximum- 
likelihood fit is performed to extract maximal informa- 
tion from the data. The fitting model accounts for in- 
terferences among different intermediate states, as well 
as the spin-dependent angular distributions of the final 
state. The large sample size, combined with the clean 
environment afforded by the presence of a J/ip or ip' in 
the final state, makes it possible to distinguish the dif- 
ferent kaon excitations that contribute to the K + it + tt~ 
final state. The results provide information not only on 
intermediate-state interactions but also on the structure 
of the kaon spectrum. By performing an additional fit in 
which the mass and width of the iiTi (1270) are floated, 
we measure the mass and width of the i^i(1270). 



1 Charge-conjugate modes are always implicit. 

2 In 2001, the Belle Collaboration measured the branching fraction 
for B+ -> J/V>-ftTi(1270) with 2% of the data presented here. 
The K+tt+tt- final state in B+ J/ipK+ 7r+7r was found 
to be dominated by the i<"i(1270), and no other structure was 
detected pj. 



Identifying the kaon excitations involved in B + — > 
J/i[iK + Tr + Tr~ and B + — > ip' K + tt + it~ can lead to a better 
understanding of the underlying theory. For example, the 
breaking of SU (3) flavor symmetry mixes the l 3 Pi and 
l 1 /*! states of the kaon system into the physical states 
if i (1270) and Ai(1400) as 

Ai(1270) = A(l 3 Pi) sin6»K + A^Pi) cos9 Kl (1) 
Ai(1400) = A(l 3 Pi) cos 6>if - K^P^wxOk, (2) 

where 9k is the 3 Pi- 1 Pi mixing angle. The value of 
9k can be related to the masses of the Ai(1270) and 
Ai(1400), to the strong decays of the Ai(1270) and 
Ai(1400), and to rates of weak decays to final states in- 
volving the Ai(1270) and Ai(1400) The measure- 
ments presented here can lead to a better determination 
of 9 K - 

The data sample used in this study was produced by 
the KEKB asymmetric-energy e + e~ collider Q and re- 
constructed by the Belle detector 0. It corresponds to 
an integrated luminosity of 492 fb _1 accumulated at the 
T(4S') resonance and contains 535 x 10 6 BB meson pairs. 



II. THE BELLE DETECTOR 

The Belle detector [fj| is a large-solid-angle magnetic 
spectrometer. A silicon vertex detector surrounds the 
interaction point and reconstructs decay vertices. A 
50-layer central drift chamber (CDC) provides charged- 
particle tracking over the laboratory polar-angle region 
17° < 9 < 150°, which corresponds to 92% of the solid 
angle in the T(4S I ) rest frame. A system of 1188 aerogel 
Cherenkov counters (ACC) and an array of 128 time-of- 
flight counters provide particle identification. An elec- 
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tromagnetic calorimeter, comprising 8736 CsI(Tl) crys- 
tals, records the energy deposited by photons, leptons, 
and hadrons. These subdetectors are surrounded by a 
superconducting solenoid 3.4 m in diameter and 4.4 m 
in length, which produces a 1.5-T magnetic field paral- 
lel to the positron beam. An iron flux return installed 
outside the coil is instrumented with large-area resistive- 
plate counters to identify muons and Kl mesons. Monte 
Carlo (MC) simulation^ are used to determine the ac- 
ceptance of the detector for the processes of interest. 

III. EVENT SELECTION 

Electron candidates are identified by combining infor- 
mation from the CDC, electromagnetic calorimeter, and 
ACC. Muon candidates are identified by extrapolating 
charged-particle tracks from the silicon vertex detector 
and CDC into the -Kl/a* detector. To identify charged 
hadrons, momentum measurements from the CDC are 
combined with velocity information from the time-of- 
flight counters, ACC, and CDC (dE/dx) @. The kaon 
identification efficiency is above 80%, while the probabil- 
ity of misidentifying a pion as a kaon is below 10%. 

Low-momentum charged tracks that curl up in the 
CDC can be reconstructed multiple times by the track 
finder. To ensure that no track is included more than 
once, criteria similar to those of Ref. [lt^ are used@ 

In studying a mode that has a J/ip or ip' in the final 
state, a key strategy is to reconstruct the J/ip only in 
its decays to e + e _ or fi + fi~, and the ip' only in its de- 
cays to e + e~, or J/ipir + ir~ . Although this choice 
abandons all but 12% of J/tp y s and 5% of ^i"s, it re- 
duces continuum backgrounds to a negligible level. The 
decays J/ip — > and ip' — > /i + /J~ are reconstructed 
by combining oppositely-charged muon candidates. The 
invariant-mass distribution is then fitted, modeling the 
J/tp and ip' as double Gaussians. Muon pairs are dis- 
carded unless they have an invariant mass within ±3cr of 
the fitted J/ip or ip' mean, where a is the width of the 
narrower Gaussian. Similarly, J/ip — > e + e~ and ip' — > 
e + e~ decays are reconstructed by combining oppositely- 
charged electron candidates. To account for energy losses 
due to final-state radiation or bremsstrahlung in the de- 
tector, any photons detected within 50 mrad of the ini- 
tial direction of an electron candidate are also included 
in the e + e~ invariant-mass calculation. Electron pairs 
are discarded unless they have an invariant mass within 
the range extending from — 4tr to +3<r of the fitted J/ip 
or ip' mean. This mass window is asymmetric about the 
mean so as to include the radiative tails of the J/ip and 



3 Lists of four-vectors for a given decay chain are generated us- 
ing EvtGen @|. The detector response is then simulated using 
GEANT §], combining randomly-triggered data with the simu- 
lated events. 

4 See Ref. [ll] for a detailed description of the event selection. 



ip', which are not completely recovered by the photon 
addition. 

Lepton track pairs that survive the mass requirements 
are fitted to a common vertex, which is constrained 
within errors to the measured interaction point. This 
vertex is then fixed, and another fit is performed, con- 
straining the dilepton invariant mass to the nominal J / ip 
or ip' mass. Since the observed widths of the J/ip and 
ip' are dominated by measurement error, this procedure 
improves the mass resolution of the B candidate. 

To reconstruct ip' — > J/ipir + iT~ decays, leptonic J/ip 
candidates are combined with a pair of oppositely- 
charged tracks that satisfy pion-identification criteria. As 
the 7T + 7r - invariant-mass distribution in ip' — » J/ipir + Tr~ 
decays is known to peak at high values [12j . the dipion 
invariant mass is required to be greater than 0.4 GeV/c 2 . 
Unless they have an invariant mass within ±3cr of the fit- 
ted ip' mean, ip' — > J/ip-K +r K~ candidates are discarded. 

To reconstruct B-meson candidates, each J/ip or ip' 
candidate is combined with a kaon candidate and two 
oppositely-charged pion candidates. Kaon and pion can- 
didates are charged tracks that satisfy identification crite- 
ria for kaons and pions, respectively, and have an impact 
parameter with respect to the fitted dilepton vertex of 
\dr\ < 0.4 cm and \dz\ < 1.5 cm i Any pion candidate 
that is identified as the product of a Kg — > tt + it~ decay 
is discarded!! 

A. B-Meson reconstruction 

Two kinematic variables can be used to identify B 
mesons. First, the reconstructed mass of a true B meson 
is likely to fall near the nominal B mass. Second, as B 
mesons are produced in the reaction 

e+e- -> T(45) -> B + B~, (3) 

the energy of each B in the T(45 l ) frame is half the to- 
tal energy of the electron and positron beams in this 
frame. Since beam-energy drifts can cause the mass of 
the T(45) to vary, it is customary to recast these kine- 
matic variables in forms that are readily corrected for 
drifts in the beam energy — namely, the energy difference 
AE and beam-constrained mass Mb c : 

AE = E*{B)~E* hc&m , (4) 

Mbc = y/Kl^-P^B)- (5) 



The impact-parameter requirement is not applied to the pions in 
ip' — > J/i/jtt + tt~ . 

To reconstruct Kg — > n + ir~ , oppositely-charged pion candidates 
are combined, and the Kg selection criteria of Ref. [0j are ap- 
plied. Both pions are vetoed if their combined invariant mass 
lies between 0.482 GeV/c 2 and 0.510 GeV/c 2 , which corresponds 
roughly to a region extending from — 4<r to +3<r around the nom- 
inal Kg mass. 
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Here, P*(B) is the momentum of the B candidate in 
the T(45) frame, while -E^ cam is half the energy of the 
T(4S f ) and is measured independently. For a correctly- 
reconstructed B meson, AE peaks at zero, and M^ c 
peaks at the nominal B mass. 

In the case of a multiparticle final state such as 
J/ipK + ir + iT~ or tp' K + tt + tt~ , multiple B candidates can 
pose a challenge. If a correctly-reconstructed B candi- 
date includes a low-momentum pion, then an additional 
B candidate can be formed by replacing that pion with 
a low-momentum pion from the other B. As the ex- 
change does not significantly affect the energy or mo- 
mentum of the B candidate, both candidates can sat- 
isfy AE and Mb c criteria. Multiple candidates can spoil 
branching-fraction measurements and distort observed 
mass spectra. To ensure that no B decay is counted 
more than once, a best-candidate selection is performed. 
First, B candidates are required to have \AE\ < 0.2 GeV 
and M bc > 5.27 GeV/c 2 . This leaves 25% of B+ -> 
J/ipK+Tt+n- events and 34% of B+ -> ip'K+n+TT- 
events with multiple candidates; these events have a 
mean multiplicity of 2.4 and 2.7, respectively. If a given 
event has multiple B candidates with the same final state, 
the charged tracks that make up each B candidate are 
fitted to a common vertex. The candidate whose ver- 
tex fit has the smallest x 2 is selected. According to MC 
studies, this procedure identifies the correct B candidate 
in approximately 55% of cases where there are multiple 
candidates. 

In the case of B + — > J /%PK + tt + tt~ , the decay B + — > 
ijj'K + is vetoed by rejecting all B candidates that have 
a J/-07r + 7r~ invariant mass between 3.675 GeV/c 2 and 
3.695 GeV/c 2 [3 According to MC studies, 1.1% of B+ -> 
ip'K + events in which the ip' decays to J/-07r + 7r~ and the 



J/tp decays to e + e" 



or /!+/! 



survive this veto. 



B. Signal and sideband regions 



The sideband normalization factor is given by 

/ pbkg dAE 

„ signal 



/ pbkg dAE ' 

sideband 



(8) 



where Pbkg is the polynomial representing the back- 
ground. The fraction of signal-region events that are 
background is estimated as 



n B 



B 
~S 



(9) 



where S and B are the numbers of events in the signal 
and sideband regions, respectively. 
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FIG. 1. Data AE distributions for £T 

and B + —>■ ip' K + tv + tv~ (bottom). The curves show the results 
of the fits described in the text. Dashed and dotted lines 
indicate the signal and sideband regions, respectively. 



Data distributions of AE for B + — > J /ipK + Tr + ir~ and 
B + — > ip' 'K + tt + tt~ are shown in Fig. [TJ The signal is 
modeled as a double Gaussian with a single mean, fixing 
the width and relative height of the wider Gaussian to 
the results of a MC fit. The background is modeled as 
a first-order polynomial. Based on these fits, the signal 
region is defined as 



30"A_E < AE - flAE < +3cta_E, 



(6) 



where fiAE is the mean of the signal peak, and <jae is 
the width of the narrower Gaussian. The sideband re- 
gion, which is used to estimate the background under 
the signal, is defined as 



0.13 GeV < AE 
0.05 GeV < AE 



Mab < -0.05 GeV, 
Hae < 0.13 GeV. 



(7) 



The small contribution from B + — ¥ X(3872)K + is not vetoed. 



IV. COORDINATE TRANSFORMATIONS 

The data in the sideband region are used to model the 
background in the signal region. Figure [5J which shows 
the distribution of M(Kmr) for B + -> J/ipK + ir + TT~ 
events in the signal and sideband regions, reveals a prob- 
lem: signal and sideband data have different end points 
in M(Ktth). 

Plotting AE versus M(Ktttt) reveals the cause of the 
discrepancy. As Fig. [3] demonstrates, the kinematically 
allowed range of M(Kirn) depends on AE. While the 
minimum value of M(Ktttt) is M(KTnr) m i n = Mk+2M^, 
the maximum value, which is attained when both the 
K + tt + tt~ system and the J/ip are at rest in the B- 
candidate's rest frame, varies with AE as M(KTnr) max = 
AE + Mb-Mj/^. Here, M K , M w , M B and Mj,^ stand 



for the nominal masses of the subscripted particles. 
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2 2.2 2.4 

M(Kjiti) (GeV/c 2 

FIG. 2. Invariant mass of the K + tt + tt~ system in B + — > 
J/ipK + iv + -K~ data. Open and filled histograms show events 
in the signal and normalized sideband regions, respectively. 
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FIG. 3. AE versus M(Ktttt) for B + ->■ J/^K+tt+tt- data. 
Dashed lines outline the signal region, and dotted lines outline 
the sidebands. The solid lines indicate the minimum and 
maximum values of M(Ktttt). 



Transforming M (Ktttt) as follows removes its depen- 
dence on AE: 

M' (Ktttt) = M(Kmr) min 

+ [M(Ktttt) - M(KTnr) min ] 
M(KTrn)° max - M{K^) mh 



M (Ktttt) 



M^TTTt), 



(10) 



Mb — Mj/^ is the value of 



Here, M^tttt)^ 

M(i4T7r7r) max at AE = 0. Figure S] shows AE versus 
the transformed coordinate M' (Ktttt). While the min- 
imum value of M(Ktttt) is unaffected by the transfor- 
mation, the maximum value is changed such that the 
maximum of M' (Ktttt) at any AE is equal to the maxi- 
mum of M(Ktttt) at AE = 0. The range of M(Ktttt) is 
compressed for positive values of AE and stretched for 
negative values of Ai?0 
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FIG. 4. AE versus M' (Ktttt) for B+ -> J/tPK+tt + tt~ data. 
The lines are defined as in Fig. [3] 



Figure [5] shows M' (Ktttt) for events in the signal and 
sideband regions. The problem of Fig.[5]has been solved: 
the end points of the transformed signal and sideband 
distributions match. 




2 2.2 2.4 
M'(Kjui) (GeV/c 2 ) 

FIG. 5. Transformed invariant mass of the K + tt + tt~ system 
111 B + J/%I)K + 7r + 7r data. Open and filled histograms 
show events in the signal and normalized sideband regions, 
respectively. 



An important feature of the transformation is that it 
does not change M(Ktttt) at AE — 0. Thus, although 
sideband and signal regions are both transformed, the 
change is minimal in the signal region. 

Just as the range of M(Ktttt) depends on AE, the 
ranges of M(Ktt) and M(tttt) also depend on AE. To 
correct for this dependence, transformations similar to 
Eq. [TU]are applied. The variable M(Ktt) is transformed 



Although correctly-reconstructed B mesons should have AE 



on average, systematic errors shift the observed mean of the 
signal AE peak away from zero by 1-2 MeV. For simplicity of 
presentation, this mean is assumed to be zero in the equations of 
this section. In fact, just as the signal and sideband regions are 
centered around the measured AE mean, the transformations 
are also made about the measured AE mean. 
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M'(Kn) = M(X7r) min 

+ [M(Kn) - M(K7r) min ] 
< M(KTrn)° max - M (K7r7r) min 
X M(X7T7r) max - M(K7T7r) min : 

and the variable M(tttt) is transformed as 

M'{im) = M(7r7r) min 

+ [M(7T7r) -M(7T7r) mln ] 

M(K7T7r)^ ax - M(i^) min 



> 

CD 

o 



(11) 



M(Kww) 



M(Kitit) y _ 



(12) 



Here, M(Kn) mill = M K + A'U, and M(irir) min = 2M V . 
Figures 111 and [3 show AE versus M(Kn) and M'(Kn). 
A similar effect is observed for M(irir) and M'(tttt). 




M(Kji) (GeV/c") 



FIG. 6. AS versus M(if7r) for B+ ->■ J/ipK + n + Tr~ data. 
The lines are defined as in Fig. [3] The concentration of 
events near 0.9 GeV/c 2 represents random particle combi- 
nations containing real _R"*(892)s. 

As Fig. [7] illustrates, transforming the M (Ktt) coor- 
dinate distorts the shapes of the if* (892) and D° back- 
grounds. This distortion must be taken into account in 
parametrizing the background for the three-dimensional 
fits of Sec. ED (i-e., in Eqs. [M] and HI . Modeling the 
distortion is straightforward. First, the peak is described 
as a Breit-Wigner or Gaussian in the untransformed co- 
ordinate, M(Kir). Using Eq. QTJ M(Kir) is then written 
as a function of M'(Kir) and AE. The expression is 
numerically integrated over the relevant range of AE to 
obtain the shape of the peak as a function of M'(Kir). 
Figure [5] demonstrates the transformation of the if* (892) 
background shape. 

The data also contain Kg and p° backgrounds, albeit 
less prominently. The distortion of these peaks by the 
M(irir) transformation is modeled by describing the K s 
as a Gaussian and the p° as a Breit-Wigner in M(7T7t), 
expressing M(tttt) as a function of M'(tttt) and AE with 
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FIG. 7. AE versus M'(Kty) for B 4 
The lines are defined as in Fig. [3] 
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FIG. 8. M(Kty) (top) and M'(Kir) (bottom) distributions of 
the K* (892) peak in generic-MC sidebands. The curve in the 
top plot is the result of a fit to a Breit-Wigner plus a polyno- 
mial background. This curve is transformed as described in 
the text and is then superimposed on the M'(Kn) distribu- 
tion in the bottom plot. 



the help of Eq. [I2l and integrating this over the appro- 
priate region of AE. 

As the main source of background in this analysis 
is misreconstructed _B-meson decays, the transforma- 
tions were checked by analyzing a generic-MC simula- 
tion of T(4S*) decays to B + B~ and B°B°, with all 
known decay modes included. The M(Ktttt), M(Ktt), 
and M(7T7r) distributions were found to display the same 
A E- dependence in MC simulation as in data. Exclud- 
ing signal events from the MC sample, the distributions 
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of background events in the signal and sideband regions 
were compared with and without the transformations. 
The transformed sidebands were found to reproduce the 
shape of the background in the signal region more ac- 
curately than the untransformed sidebands, especially at 
high M(Kmr), M(Kn), and Af(7r7r)H 

Although some discrepancy was observed between the 
background in the signal and sideband regions near the 
Kg and p masses, this is mostly independent of the trans- 
formation and is taken into consideration when calculat- 
ing systematic errors. 

The transformations of Eqs. [T0llT2l were also applied to 
B + -> ip'K+TT+n- , with M J/M , replaced with M^> . The 
results of the checks were the same. 

For simplicity, the variables M'(Ktttt), M'(Ktt), and 
M'(7T7r) are henceforth referred to as M(Ktttt), M(Ktt), 
and M(7T7r), respectively. 



V. TOTAL BRANCHING FRACTIONS 

Branching fractions for S-meson decays to 
J/ipK + ir + ir~ and ip' K + tt + it~ final states are mea- 
sured using a background-subtraction technique^ For 
each final state, data events in the signal and sideband 
regions are distributed into cubic bins in M 2 (KTrir), 
M 2 (Kn), and M 2 (irir). The number of signal events 
observed in each bin is calculated as 

Ni = Si-f B Bi, (13) 

where Si and Bi are the numbers of signal-region and 
sideband-region events, respectively, that fall into the ith 
bin, and fs is the sideband normalization factor given by 
Eq. [SJ The fraction of charged B mesons that decay to 
the final state in question can be expressed as 

where Si is the signal efficiency in bin i, and N B is the 
total number of charged B mesons in the data sample. 
Assuming equal rates for T(45) -)• B + B~ and T(4S) -> 
B°B°, Nb is equal to the number of B pairs produced, 
which is measured independently. 

To determine the signal efficiency, we generate 10.7 x 
10 6 nonresonant B ± decays to each of the two final states 
of interest. We then reconstruct these signal-MC events, 
applying the same event-selection requirements as with 
data. We bin the generated events according to the gen- 
erated values of M 2 (KTnr), M 2 (Ktt), and M 2 (ttit), and 
the reconstructed events according to the reconstructed 



values of M 2 (Kirir), M 2 (Ktt), and M 2 (iytt). The effi- 
ciency in each bin is the ratio of reconstructed to gen- 
erated events in that bin. Figure |9] shows the depen- 
dence of the efficiency on the three variables. Figure [TU1 
shows the corresponding data distributions. The overall 
efficiency is (19.85 ± 0.01)9^3 for B+ -> J/^K+tt+tt- 
and (6.58 ± 0.01)% for B+ -> tp' K+tt+tt- . The number 
of efficiency-corrected signal events observed is (4.14 ± 
0.06) xlO 4 for B+ -> J/ipK+n+n- and (1.12±0.05) x 10 4 

for B+ ->■ Ip'K+TT+TT-. 

This method of measuring branching fractions auto- 
matically corrects for efficiency variations over the phase 
space. It also makes no assumptions as to the shape of 
the signal in AE. 



A. Systematic errors 

The systematic error in the branching fractions is es- 
timated by adding in quadrature various contributions, 
which are assumed to be uncorrelated. Where possible, 
a correction is applied. 

Since we use MC simulation to determine the signal ef- 
ficiency in Eq. I14| any discrepancy in signal-reconstruc- 
tion efficiency between data and simulation will result 
in a systematic error. Based on studies of the track- 
reconstruction efficiency, we include a systematic error of 
1.0% for each lepton track, 1.4% for each pion track, and 
1.2% for each kaon track, adding linearly. Based on stud- 
ies of the lepton-identification efficiency, which show that 
the simulation underestimates the lepton-identification 
efficiency, we apply a correction factor of 0.984±0.019 for 
each electron track and 0.962±0.031 for each muon track. 
Based on studies of the kaon identification efficiency, we 
include a systematic error of 1% for B + — > J/ipK + ir + ir~ 
and 2% for B+ -> tfj' K+tt+tt" . 

The bin size of 0.15 GeV 2 /c 4 is chosen based on the 
dependence of the efficiency-corrected signal yield on the 
bin size. The error associated with this choice is taken 
to be the rms of the signal yield in the region between 
0.1 GeV 2 /c 4 and 0.2 GeV 2 /c 4 . 

In the AE distributions for B + — > J/i(>K + ir + Tr~ and 
B + — > ip' K + ir + ir~ nonresonant MC simulation, shown 
in Fig. 111! a small polynomial background can be seen 
under the peak. Since all the events in the MC sample 
include a signal decay, this "background" is made up of 
misreconstructed signal events. Although these events 
are included as signal in the efficiency calculation, they 
are removed by the background-subtraction procedure. 
The fraction of signal that is subtracted in this way is 
found to be (3.75 ± 0.91)% for B+ -» J/^K+n+n-, 
and (5.7 ± 1.2)% for B+ -)• ip'R+n+n-. The observed 



9 For details, see Ref. lilt 11 Throughout this paper, when a single error is presented, it is 

10 Peaking backgrounds are not expected in these final states and statistical; when two errors are presented, the first is statistical, 

were not seen in generic-MC simulation. and the second is systematic. 
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FIG. 9. Dependence of the signal efficiency on the kinematic variables for B + — > J /ipK + -r + -k (left) and B + — > ip'K + ir~* 
(right). 
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branching fractions are corrected for this effect, and the 
associated uncertainty is included as a systematic error. 
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the ip' — > J/ipir + ir veto in the reconstruction of B + — > 
J/ipK + Tr + Tr~ . This branching fraction is 

B(B + -> tP'K+) = (6.65 ± 0.17 ± 0.55) x 10~ 4 , 

which is consistent with the previously-measured value 
of (6.48 ±0.35) x 10" 4 QJ]. 

Our B + — > J/ipK + ir + ir~ branching-fraction measure- 
ment represents a significant improvement over previ- 
ous measurements (Tj. It is consistent with Ref. [IH 
but inconsistent with Ref. [l6| at the 3.4-ct level. Our 
B + — > ip' K + Tt + Tt~ branching-fraction measurement is 
also a significant improvement over the previous mea- 
surement uM- 



VI. AMPLITUDE ANALYSES 

To study the resonant structure of the K + Tr + Tt~ final 
state in B + -> J/tfjK + n + TT^ and B + -> tp' K + n + Tr~ , 
we perform amplitude analyses. Using an unbinned 
maximum-likelihood method, we simultaneously fit the 
data in the three dimensions M 2 (Kirir), M 2 (Kir), and 
M 2 (tttt). 



To determine the sideband normalization factor fg in 
Eq. [T31 the data AE distribution is fitted as described in 
Sec. MI Bl In this fit, the background under the signal 
is parametrized as a first-order polynomial. To estimate 
the error introduced by this assumption, a second fit is 
performed, parametrizing the background as a second- 
order polynomial. The fractional change in the signal 
yield is taken as a systematic error. 

Since the signal and sideband regions are defined based 
on the results of fitting the data AE distribution, (iae 
and & ae in Eqs.[6]and[7]are varied within the fit errors. 

In the MC sample used for determining the efficiency, 
J/ip'a from the signal B are forced to decay to e + e~ 
or /i + /Lt _ , and ip n s from the signal B are forced to 
decay to e + e~, or J / ipir + ir~~ . Thus, to ob- 

tain branching fractions for B + J /ipK + tt + ti~ and 
B + — > ip' K + tt + tt~ , the branching fractions measured us- 
ing Eq.[l4]are divided by previously- measured values [l4| 
of these J/ip and i/j' decay rates. The uncertainties of 
these previous measurements are included as a system- 
atic error. 

Finally, the error in Nb is 1.3%. Table|I]lists the com- 
ponents of the systematic errors. 



B. Results 

The measured branching fractions are 

B(B + -> J/ipK + TT + ir-) = (7.16 ± 0.10 ± 0.60) x 10~ 4 , 
B(B + -> iP'K+tt+tt-) = (4.31 ± 0.20 ± 0.50) x 10~ 4 . 

As a cross-check, we also measure a branching fraction 
for B + — > tp'K + , using a similar method but reversing 



A. Fitting technique 

Signal-region data are fitted by maximizing^ the log- 
likelihood function, which is given by 



(15) 



where the sum is over the events in the signal region, 
xi is the vector of coordinates for a given event (i.e., 
x = [M 2 (K7nr),M 2 (Kn),M 2 (7nr)}), a is the vector of 
parameters with respect to which I is maximized, and p 
is the probability-density function (PDF) that is used to 
model the observed distribution. 

The distribution of events in the signal region is mod- 
eled as 



p(x; a) = n B 



Pb(x) 



JpB(x)d 3 x ' S Jps(x;a)d 3 x 



Ps(x;a) 



(16) 



where pg and p$ describe the observed shapes of the 
background and signal, respectively. The constants ng 
and ns are the background and signal fractions in the 
signal region; the former is given by Eq. [31 and the latter 
isl-7i B E3 

The observed signal distribution ps is expressed as 



p s (x; a) = e{x)cj)(x)s{x\ a), 



(17) 



12 Standalone MINUIT [3 is used for all maximizations in this 
section. 

13 The background fraction ns is corrected for the oversubtraction 
effect described in Sec. IV Al 
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TABLE I. Components of the systematic error in the branching-fraction measurements, expressed as a percentage of the 
branching fraction. 



Component 


B + -> J/lpK + TT + TT- 


B+ ->■ tp'K+n+TT- 


B+ ->■ ip'K+ 


MC statistics 


0.18 


0.19 


1.12 


Tracking efficiency 


6.0 


8.7 


6.0 


Lepton-ID efficiency 


5.1 


5.1 


5.0 


Kaon-ID efficiency 


1.0 


2.0 


1.0 


Binning 


0.13 


0.19 


0.26 


Oversubtraction 


0.95 


1.3 


0.95 


Background shape 


1.6 


1.3 


0.19 


Signal/Sideband regions 

J /if) or branching fraction 


0.35 


4.4 


0.27 


0.71 


1.9 


0.67 


N B 


1.3 


1.3 


1.3 



where e is the detector efficiency, 4> is the phase-space 
density, and s is the raw signal function. 

Using nonresonant MC simulation, we have measured 
the detector resolution to be approximately 3-4 MeV/c 2 
in each of the three coordinates M(Kirir), M(Ktt), and 
M(7T7r). Since this is smaller than the width of any reso- 
nance included in the fits, we neglect the effect of detector 
resolution on line shapes. 

The following five sections describe the methods fol- 
lowed in performing the integrals of Eq. [16] and in ob- 
taining the functions ps(x), e(x), 4>(x), and s(x;a) in 
Eqs.[IUand[T7] 



B. Normalization procedure 

The integrations of Eq. [16] are performed numerically, 
using Simpson's rule. A step size of 0.010 GeV 2 /c 4 for 
B+ J/iJjK+it+it- and 0.005 GeV 2 /c 4 for B+ ->• 
tp' K + tt + tt~ is used in each dimension The three- 
dimensional region of integration can be determined 
by noting that the minimum and maximum values of 
M 2 (Kitit) are given by 



M 2 (Kirir) 
M 2 (Kttit), 



(M K 



- 2M„y 



(18) 
(19) 



where Mb, Mk, M n , and M^, are the nominal values of 
the subscripted particles. For a given value of M 2 (Kim), 
the minimum and maximum values of M 2 (irir) are 



M 2 (tttt\ 



(2M W ) 2 , 



M 2 (^)max = (^M 2 (KtTTt) - M K )' 



(20) 
(21) 



For given M 2 (Ktttt) and AI 2 (tttt), the minimum and 
maximum values of M 2 (Ktx) are 



M 2 (Kn)™* = - 



M 2 (Ktttt) + M 2 K + 2M 2 - M 2 (tttt) 



± 



V /l-(2M 7r /M( 7 T7r)) 5 



x ^/M 2 (Kitit) - (M K + M(nir)Y 



x y/M 2 (Kirir) - (M K - M^tt)) 2 



(22) 



Figure [T3] shows the calculated kinematic boundaries for 
B + — > J / %jj K + ix + tt~ , along with the observed distribu- 
tions of sideband data, for a slice in M 2 (Kirir). 
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a 3.5- 



2- 
1.5- 
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1^(7171) (GeV 2 /c 4 ) 

FIG. 12. Scatterplot of M 2 (Ktt) versus M 2 (im) for 
M 2 (Kitk) between 4.0 GeV 2 /c 4 and 4.5 GeV 2 /c 4 in sideband 
data for B + —t J/^K + it + it~ . Blue and red curves show the 
calculated boundaries corresponding to the low and high edge, 
respectively, of the plotted M 2 (Kitty) region. 



14 The larger step size is necessary for B~*~ — ► J/ipK + tt + tt~ because 
of the larger phase space, which significantly increases the CPU 
time required for the integration. 



Events that do not fall within the calculated bound- 
aries are excluded from the fits. Because of the coordi- 
nate transformations of Sec. IIV1 such events are rare: 36 
of the 12913 sideband events and 3 of the 10594 signal- 
region events for B + — > J/i/;K + tt + tt~ , and 3 of the 2230 
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sideband events and none of the 1176 signal-region events 
for B + — > ip' K + ir + ir~ fall outside the boundaries. 



C. Background functions 

To determine the three-dimensional shape of the back- 
ground in the signal region (i.e., pb(x) in Eq. I16[) . an 
unbinned maximum-likelihood fit is performed on the 
sideband-region data. The log-likelihood function to 
maximize is given in this case by 



(-b{o.b) = In 



PB{Xj\a.B) 

J pb(x; a B ) d 3 x' 



(23) 



where the sum is over the events in the sideband region. 
The maximization is performed by varying the parame- 
ters ob, which are then fixed at their optimal values in 
fitting the signal region. 

The background is modeled as the sum of a combina- 
torial term and a set of noninterfcring resonances. For 



B- 



J/tPK+it+tt-, 



p B (x;a B ) 



^2a xi Ti(x) 

i=0 

[o-k* (892) BW K , i892) {y) + a p BW p (z) 



Lj=o 



2 

E 

L/c=0 



a z kT k (z) 



and for B + -» 

Pb(x;S b 



+ a D Go (y) + a Ks Gr s (z) 

if)' K +r K + TT~ , 



(24) 



-2a- 



^a xi Ti(x) 

i=0 

a K ' (892) BW K , ( 89 2) (y) + a p BW p (z) . (25) 



In Eqs. [24] and [25j T n represents an nth-order Cheby- 
shev polynomial. The variables x, y, and z stand for 
M 2 (Ktttt), M 2 (Kn), and M 2 (tttt), respectively, and are 
defined over the intervals 



(M K + 2M„Y 
(M K + M n ) 2 , 
(2M„) 2 , 



^rnax 


= (M B 


£/max 


= (M B 


•^m ax 


= (M B 



M^p) 2 , 



(M B - M^ - MkY 



The peak functions 5^^.(892), Go, Gk S i an d BW P are 
obtained as described in Sec. IIVI Each peak function 
P(x) is normalized over the kinematically-allowed phase 
space to satisfy 



(26) 



1 = / e~ 2x P(x)d 6 x. 



The factor of e 2x that modulates the peak functions was 
found empirically to produce a good fit to the sideband 
dataEl 

Table |ll] lists the fitted parameters of the background 
functions. The statistical error in each parameter is de- 
fined as the change in that parameter required to reduce 
the log likelihood by 1/2. The fitted functions, normal- 
ized to the total number of events in the fit, are shown 
projected onto the three axes along with the sideband 
data in Fig. [H Figures [H and HU show M 2 (Kir) and 
M 2 (tttt) projections for slices in M 2 (Kirir). 



TABLE II. 


Fitted values of the backj 


ground-function param 


eters (Eqs. 


MlandES)). 




Parameter 


B+ -> J/tPK+iy+tt- 


B+ iIj'K+tt+ti- 


CLxO 


1.0 (fixed) 


1.0 (fixed) 


Clxl 


-1.5901 ± 0.0048 


-1.238 ± 0.036 


dx2 


0.9423 ± 0.0086 


0.480 ± 0.045 


dx3 


-0.4737 ± 0.0087 


-0.221 ± 0.021 


CLxi 


0.1778 ± 0.0067 




dx5 


-0.0488 ± 0.0033 




a y o 


1.0 (fixed) 




dyl 


0.088 ± 0.021 




a z o 


1.0 (fixed) 






-0.022 ± 0.022 






0.129 ± 0.018 




IK* (892) 


0.0353 ± 0.0065 


0.0161 ± 0.0062 


an 


0.0007 ± 0.0011 




a Ks 


0.0061 ± 0.0023 




a p 


0.086 ± 0.012 


0.0352 ± 0.0099 



As a measure of goodness of fit, a x 2 variable is cal- 
culated by distributing the data into cubic bins that are 
0.1 GcV/c 2 wide on each side. The normalized PDF, 
with the parameters set to their best-fit values, is inte- 
grated over each bin and multiplied by the total number 
of events in the fit to determine the number of events 
expected in the bin. Adjacent bins are combined until 
each bin has at least 6 data events. A \ 2 variable for the 
multinomial distribution is then calculated as 1 1 91 



X 2 = 2 



Whins 

E 



m In I — 

Pi 



(27) 



where A^ins is the total number of bins used, is the 
number of observed events in a given bin, and Pi is the 
number expected in that bin based on the PDF. 

If the expected distribution pi were obtained by a 
binned maximum-likelihood fit of the data distribution 
Hj, the number of degrees of freedom associated with this 
X 2 would be reduced by the number of fit parameters 



15 Since there are more low-energy particles than high-energy 
particles, the background peaks are more pronounced at low 
M 2 (Knir). Combining a if* (892) with a random pion, for ex- 
ample, will tend to produce a low value for M 2 (Knn). 
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FIG. 13. Results of sideband fits for B + ->■ J/^K+n+n' (left) and B + -> ^'K + n + n' (right). Data (points) and fits 
(histograms) are shown projected onto the three axes. The red histograms show the overall background functions. The 
combinatorial components are shown in gray, while the .ft'* (892), p, Ks and D backgrounds are shown in blue, magenta, cyan, 
and green, respectively. The A"" (892) and p peaks are broader in B + —¥ ip' K + -k + tt~ than in B + — > J/ipK + -K + Tv~ because the 
distortion shown in Fig. [8] is larger in the former mode. 



iVp ar and would be given by Nuof = Nbins - N pax - 1. 
If, on the other hand, the two distributions were not cor- 
related by a fit, the number of degrees of freedom would 
be iVrjoF = A/bins — !• Since, in this case, the distri- 
butions are related by an unbinned maximum-likelihood 
fit, the true ./Vdof can be expected to lie between these 
extremes (2fjj . 

For the B + — > J/ipK + Tr + Tr" sideband-data fit, £b — 
-21484.6, while X 2 = 1709.5 with N hins = 1707 and 
iV par = 12. For the B + -> ^'K+tt+tt^ sideband-data 
fit, l B = 822.7, while X 2 = 286.8 with iVbins = 294 and 

^par = 5 - 

D. Efficiency functions 

The dependence of the detector efficiency on the kine- 
matic variables (i.e., e(x) in Eq. rTTl) is obtained for three- 
dimensional bins, 0.15 GeV 2 /c 4 -wide on each side, using 
nonresonant signal-MC simulation as described in Sec.lVl 
and illustrated in Fig. |H1 The function is implemented 
as a lookup table: the efficiency for a given data point is 
the efficiency in the corresponding bin. 



E. Phase-space densities 

Four-body phase-space densities (i.e., cj>(x) in Eq. [T7|) 
for B + -t J/tpK + TT + ir~ and B + i/j'K+tt+tt^ are ob- 
tained by using GENBOD 21] to generate final-state- 
particle four-momenta that are weighted by the density 
of states in phase space [22[ . For each decay mode, 
10 8 events are generated. Event phase-space weights 
are distributed into cubic bins in M 2 (Ktttt), M 2 (Ktt), 
and M 2 (tttt), with a bin width of 0.02 GeV 2 /c 4 . The 
phase-space density is implemented as a lookup table: 
the value of <j>(x) for a given data point is the total phase- 
space weight in the corresponding bin[j£| In Fig. [TBI 
the three-dimensional histogram of phase-space weights 
is projected onto the three axes, showing the distribution 
that signal events would have in the absence of resonant 
effects. 

Figure fTTH does not indicate the functional form of 4>(x), 
since the projection onto a single dimension effectively in- 



Boundary effects are insignificant. 



13 



0.60 GeV 2 /c 4 < M (Ktc7i) < 1 .46 GeV 2 /c 4 

<J; 120- 

CD 

O 100 
o 

S 80 

° 60- 
co : 

g 40- 

lS 20 




2.5 



3 3.5 4 
M 2 (Kti) (GeV 2 /c 4 ) 



0.60 GeV 2 /c 4 < M z (Kjiti) < 1 .46 GeV z /c 4 

Q 




1.5 



2 2.5 

M 2 (roi) (GeV 2 /c 4 ) 




2.5 



^ 140 
> 

CD 120 
O 

100 
■* 

9 80 
o 

» 60 

1 40- 





3 3.5 4 
M 2 (Kti) (GeV 2 /c 4 ) 





1.5 



2.32 GeV 2 /c 4 < M 2 (KitJt) < 3.18 GeV 2 /c 4 
"o 200 



£ 180i I ■ \\ I 

O 160- 'inj 1 , 



2.5 



M 2 (3i3i) (GeV 2 /c 4 ) 




M 2 (Kti) (GeVVc 4 ) 



M 2 (3i7i) (GeV 2 /c 4 ) 



100 




M 2 (Kti) (GeVVc 4 ) 



4.04 GeV 2 /c 4 < M 2 (K;t7i) < 4.90 GeV 2 /c 4 
50- 




3 3.5 4 
M 2 (Kti) (GeV 2 /c 4 ) 




M 2 (3i3i) (GeV 2 /c 4 ) 



4.04 GeV 2 /c 4 < M 2 (Kjiti) < 4.90 GeV7c 4 



CD 

O 
o 
o 

O 

CO 
CD 




2 2.5 

M 2 (roi) (GeV 2 /c 4 ) 



FIG. 14. B + — > J/tpK + TY + Tr sideband data (points) and fit results (histograms) for slices in M 2 (Kirn) . The fit components 
are color coded as in Fig. 1131 
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FIG. 15. B + — > ip'K + iT + ir sideband data (points) and fit results (histograms) for slices in M 2 (Kttit). The fit components 
are color coded as in Fig. 1131 
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tegrates over the other two dimensions, and the region of 
integration is the complicated one described in Sec. I VI Bl 
In Fig. [17l the same projections are performed over a nar- 
row slice in each of the other two dimensions, to illustrate 
the dependence of the function <f>(x) on each variable. 



F. Signal functions 

The K + tt + tt~ final state is modeled as a nonresonant 
signal plus a superposition of initial-state resonances R\ . 
The latter are assumed to decay through intermediate- 
state resonances i? 2 as i?i — > ai? 2 , Ri — > be, where a, 
6, and c arc the final-state particles. Specifically, the 
function s(x; a) of Eq. [T7]is expressed as 



s(x;a) = s(x;a k ) 

— |ci nr yl nr (x) | -\- ^ 
Ji 



^2aj lJ2 Aj lJ2 (x) 



(28) 



Here, J\ and J 2 stand for the spin-parity (J p ) of Ri and 
i? 2 , respectively. Resonances with different Ji are added 
incoherently, while those with the same J± are added co- 
herently. The parameters varied in the fit are the com- 
plex coefficients a nr and aj ± j 21 collectively referred to as 
dfe. While the nonresonant signal is assumed to be con- 
stant over the phase space, 



A nr {x) — 1: 



(29) 



the resonant decay amplitudes Aj 1 j 2 are expressed as 

A JlJ2^) = a JlJ2 



l abc 



iM Rl T Rl 



y/M R2 T R2 (mb c ) 



Ml 

K2 



mt 



%c ~ iM R2 T R2 {m bc ) 
where T R2 (m bc ) is the mass-dependent width 

2J 2 + 1 / n,r 

I M R2 



r_R 2 (mf, c ) = T R2 ( — 



F, 



and F R is the Blatt-Weisskopf barrier factor 



1 



Vi + FV 



for J 2 
for J 2 



R- 

0, 
1, 



(30) 



(31) 



(32) 



V9 + 3i? 2 g 2 + i^4 
y/9 + 3R 2 q 2 + i?V 



for J 2 = 2. 



The meson radial parameter R is set to 1.5 (GeV/c) -1 . 
The function aj x j 2 describes the spin-dependent angu- 
lar distribution of the final state and is shown for vari- 
ous combinations of J\ and J 2 in Table IIIII Resonances 
with spin greater than two are not included in the fitting 



TABLE III. Angular distribution of the K + ir + ir~ final state 
for various combinations of initial and intermediate-state spin 
parities. See Ref. |23| for derivation and conditions of applica- 



bility. 


Ji 


J2 


U J1 J2 


Any 


0+ 


1 


0+ 

0- 


+ 

T— 1 1— 1 


(l + z 2 )cos 2 <9 




1" 
1 + 


1 + z 2 cos 2 6 




1 + 

1" 


1 - cos 2 6 


1 + 


2+ 




1" 


2" 


X li + ocos U + \)Z (COS u — 1/oJ J 


2+ 

2- 


+ 1 


3 + (1 + 4z 2 ) cos 2 6 


2+ 


1" 


1 - cos 2 6> 


2~ 


1 + 




2+ 


2+ 


l + z 2 /9 + (2 2 /3-l)cos 2 


2" 


2" 


- z 2 (cos 2 6» - 1/3) 2 


2+ 
2- 


2~ 
2+ 


1 + z 2 /3 + z 2 cos 2 6 + z 4 (cos 2 e ~ 1/3) 2 



model. In cases where there is more than one covariant 
spin amplitude, only the lowest spin is included. 

In Ecis. l3"01l3"2"l and Table Unl the nominal masses of the 
resonances R\ and i? 2 are denoted by M Rl and M R2 , 
and the nominal widths by T Rl and T R2 . The angle 9 is 
between a and b in the be rest frame and can be expressed 
as 



cos 9 



m bc 



4:pqm abc 



( m lbc ~ ™ 2 a )(m 2 - m 2 ) 



m 



be 



The variable z is given by 

z = p/m abc . 



(33) 



(34) 



The breakup momentum p is the momentum of a or be 
in the abc rest frame: 

2 _ ( m lbc - ( m a + m bc) 2 )(m 2 abc - (m a - m bc ) 2 ) 



p = 



Am 2 



while q is the momentum of b or c in the be rest frame: 



(m 2 c - (m b + m c ) 2 )(m 2 c - (m b - m c ) 2 ) 



4m 2 c 



(36) 



where the constant go is the value of q evaluated at m bc — 
M R2 . 

Since the components of the signal function are not 
individually normalized, it is not meaningful to compare 
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FIG. 16. Projections of the three-dimensional histogram of phase-space weights onto the three axes for B + — s- J/ipK + ir + -K 
(left) and B+ -> ^'K+tt+tt- (right). 
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FIG. 17. Phase-space densities for selected regions, as a function of M 2 (Kttit), M 2 (Ktv), and M 2 (7T7t), for B + — > J/tpK + TT + TY 
(left) and B + — > -i/j'^+tt+tt - (right). In each case, the region selected is indicated above the plot. 
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the moduli of the complex coefficients a k in Eq. [28l A 
decay fraction is therefore calculated for each component 
by integrating the component over the kinematically- 
allowed region and dividing by the integral of the full 
signal function 

= f <b(x)\a k A k (x)\ 2 d 3 x 
J <fi(x)s(x; a)d 3 x 

The integrations in Eq. [37] are performed as described in 
Sec. I VI Bl Because of interference effects, decay fractions 
for a given final state will not, in general, add up to unity. 

G. Statistical errors 

As with the sideband-region fits, the statistical uncer- 
tainties in the fit parameters (i.e., moduli and phases) 
are determined by the fitter: the error in a given param- 
eter is the change in that parameter that reduces the log 
likelihood by 1/2. The statistical uncertainties in the de- 
cay fractions, on the other hand, are more complicated. 
Since a given decay fraction involves the integral of the 
full signal function, the error in a single decay fraction 
incorporates the errors in all of the parameters. To de- 
termine the statistical errors in the decay fractions, 1000 
sets of correlated signal-function parameters are drawn 
from Gaussian distributions using the fitted parameter 
values and the error matrix^ Decay fractions are cal- 
culated for each set of generated parameters. The rms 
of the resulting distribution provides an estimate of the 
statistical error in the decay fraction. 

H. Systematic errors 

Several sources of systematic error are considered, as 
described below. They are added in quadrature to obtain 
the systematic errors reported in Sec. I VI II 

1. Background parametnzation 

A possible source of systematic error in the fits is the 
fixed background fraction in Eq. [TH While the er- 
ror in riB is small, the correction for the oversubtraction, 
described in Sec. IV Al lowers tib by 10.8% for B + —> 
J/iiK+iT+n- and by 11.4% for B+ -> ip'K+ir+n- . The 
systematic error associated with this correction is esti- 
mated conservatively as the change in each parameter 
when the fits are performed with the uncorrected values 
of ub- 

There may be an additional systematic error if 
the background in the signal region is not correctly 



parametrized by the shape determined by fitting the side- 
bands. As noted in Sec. IIV| generic-MC studies suggest 
that not enough of the Kg and p background peaks are 
removed by the sideband subtraction. To estimate this 
error, a fit is performed in which the coefficients of the 
background peaks in Eqs. l24l and l25l are doubled. 



2. Efficiency 

To estimate the error introduced by binning the effi- 
ciency information, the fits are repeated using bin sizes 
of 0.10 GeV 2 /c 4 and 0.20 GeV 2 /c 4 for the efficiency. The 
average absolute change in each parameter is the estimate 
of the error. 

Another possible source of error is that the MC sim- 
ulation may not faithfully reproduce the detector effi- 
ciency for low-momentum particles. To test for such an 
effect, two additional fits are performed. In the first fit, 
only charged particles with a momentum greater than 
200 MeV/c are included. In the second fit, the \dr\ and 
\dz\ requirements described in Sec. IIIII are loosened from 
0.4 cm to 0.8 cm, and from 1.5 cm to 3.0 cm, respec- 
tively. The changes in each parameter observed in these 
two fits are added in quadrature to obtain an estimate of 
the error due to inaccuracies in the efficiency estimation. 

Using only the three variables M 2 (Kirir), M 2 (Kw), 
and M 2 (mr) in this analysis is equivalent to integrating 
over variables that describe the relative momentum of the 
J /if) or ip' with respect to the K + n + ir~ system. In this 
integration, the terms corresponding to K + it + tc~ states 
with different initial-state spin-parity cancel out, produc- 
ing Eq. [23 This cancellation, however, is exact only if 
the detector efficiency is flat over the extra variables. To 
determine the effect of neglecting these variables, an ad- 
ditional set of fits is performed, in which the efficiency 
in Eq. [TjJ is calculated as a function of the two angles 
between the J/ip or -0' and the K + ir + n~ system, rather 
than M 2 (Knn), M 2 (Kn), and M 2 (mr). The resulting 
fitted parameters are compared to those obtained by a 
fit in which the efficiency is held constant The abso- 
lute change in each parameter is found to be small (less 
than 15% of the statistical error) and is included in the 
systematic error. 



3. Integration step size 

To estimate the error introduced by the finite step size 
used in the numerical integrals of Sees. IIVI and IVIBI 
the fits are repeated, using a step size of 0.005 GeV 2 /c 4 
for B+ -> J/i/jK+ir+ir- and 0.010 GeV 2 /c 4 for B+ -> 



Correlated Gaussian distributions are generated using CORSET 
and CORGEN [H). 



If the efficiency is calculated as a function of all five dimensions, 
the accuracy of the result becomes dominated by the MC statis- 
tics. 



18 



TABLE IV. Masses, widths, and spin-parity values of the reso- 
nances included in the fits. With the exception of the 7^(1600) 
parameters (discussed in Sec. I VI J"3)) . all values are from 



IXcoUIldjIlCfc; 


Mass 
(MeV/c 2 ) 


Width 

(MeV/c 2 ) 


I p 


P° 


775.49 


± 


0.34 


146.2 


± 


0.7 


1" 


u> 


782.65 


± 


0.12 


8.49 


± 


0.08 


1" 


/o(980)° 


980 


± 


10 


50 


± 


■iii 

10 


0+ 


/a(1270)° 


1275.1 


± 


1.2 


185.1 


± 


2.9 
2.4 


0+ 


K* (892)° 


896.00 


± 


0.25 


50.3 


± 


0.6 


1" 


A"i(1270)+ 


1272 


± 


7 


90 


± 


20 


1+ 


A"i(1400)+ 


1403 


± 


7 


174 


± 


13 


1+ 


K*(UW)+ 


1414 


± 


15 


232 


± 


21 


1" 


^o*(1430)+ 


1425 


± 


50 


270 


± 


80 


0+ 


if 2 *(1430)+ 


1425.6 


± 


1.5 


98.5 


± 


2.7 


2+ 


if 2 *(1430)° 


1432.4 


± 


1.3 


109 


± 


5 


2+ 


#(1600)+ 


1605 


± 


15 


115 


± 


15 


2" 


iT(1680) + 


1717 


± 


27 


322 


± 


110 


1" 


A" 2 (1770)+ 


1773 


± 


8 


186 


± 


14 


2" 


if 2 *(1980)+ 


1973 


± 


26 


373 


± 


68 


2+ 



ip' K + ir + n~ . The change in each parameter is an esti- 
mate of the uncertainty associated with the numerical 
integration. 
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4- Modeling of the signal 

The masses and widths of the resonances included in 
the fits are listed in Table IIVI To estimate the system- 
atic error associated with the uncertainties in these quan- 
tities, the fits are repeated, varying each fixed quantity 
within its errors. For each mass or width, the average 
absolute change in each parameter is recorded. These 
average changes are then added in quadrature. 

In fitting the B + — > J/ipK + ir + ir~ data, the modulus 
for iq(1430) K*(892)tt is allowed to float. Relative 
to this modulus, the moduli for if|(1430) -> Kp and 
if! (1430) — > Kuj are fixed based on previously-measured 
relative branching fractions [14]. To estimate the associ- 
ated systematic error, additional fits are performed, vary- 
ing these branching fractions within their uncertainties. 

I. Results 

Table [V] lists the values of the moduli and phases of 
the complex coefficients a,k of Eq. [28] obtained by fitting 
signal-region data for B + — > J/ijjK + ir + ir~ , as well as 
the corresponding values of the decay fractions, given 



The phases of the three submodes are allowed to float. 



FIG. 18. Results of signal-region fits for B + — > J/ip K + -k + -k~ . 
Data (points) and fits (histograms) are shown projected onto 
the three axes. The fit components are color coded as shown 
in Fig. [H 



by Eq. [37] The fitted PDF is shown projected onto 
the three axes, along with the data, in Fig. [T5] Fig- 
ure [22] shows M 2 (Kir) and M 2 (ttit) projections for slices 
in M 2 (Kirir). The legend is presented in Fig. Ql|l In this 
fit, i = -10575.4, while x 2 = 1475.9 with iVbins = 1202 
and iV par = 28. 

Similarly, Table I VI I shows the fitted parameters for 
B + — > tp' K + ir + ir~ signal-region data, as well as the 
corresponding decay fractions. Figure [5D] shows the fit- 
ted PDF and data projected onto the three axes, while 
Fig. |2"B1 shows M 2 (Kir) and M 2 (irir) projections for slices 
in M 2 (Ktttt). In this fit, t = 638.3, while X 2 = 180.1 with 
N hins = 168 and iV par = 10. 

Finally, the B + — > J/ipK + ir + iT~ signal-region data are 
fitted again, this time floating the mass and width of the 
.Ki(1270). The fitted mass and width are 

M Kl (i2io) = (1248.1 ± 3.3 ± 1.4) MeV/c 2 , (38) 
IV^mo) = (H9.5 ± 5.2 ± 6.7) MeV/c 2 . (39) 

Table IVIII shows the fitted parameters, along with the 
corresponding decay fractions. Figure [5T] shows the fit- 
ted PDF and data projected onto the three axes, while 
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FIG. 19. Legend for Figs. [TS] [23 HOI OH HU and [Ml 
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FIG. 20. Results of signal-region fits for B + — > ip' K + iy + tt~ . 
Data (points) and fits (histograms) are shown projected onto 
the three axes. The fit components are color coded as shown 
in Fig. [191 



Fig.HMlshows M 2 (Kir) and M 2 (tttt) projections for slices 
in M 2 (Ktttt). In this fit, t = -10525.3, while \ 2 = 
1404.4 with TVbins = 1202 and N pal = 30. 

A comparison of TablesMand lVlTl reveals that in many 
cases, the effect of floating the mass and width of the 
iiTi (1270) results in a substantial decrease of the sys- 
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FIG. 21. Results of signal-region fits for B + —> ,J/ipK + -K + ir~ , 
with the mass and width of the JsTi(1270) floated. Data 
(points) and fits (histograms) are shown projected onto the 
three axes. The fit components are color coded as shown in 

Fig. [[I 



tematic error, which is somewhat offset by an increase 
in the corresponding statistical error. In particular, the 
fCi(1270) — » K*(892)ir decay fraction is especially sensi- 
tive to the JsTi (1270) mass and width. 

In any fit involving many floating parameters, local 
likelihood maxima can present a problem. To ensure 
that the fit results are global maxima, 100 additional 
fits were performed for each of the three cases, select- 
ing random starting values for the parameters. None of 
these fits yielded better likelihoods than those presented 
above. The local maxima encountered in the course of 
this test are discussed in the Appendix. 



J. Discussion 

1. Signal components 

In choosing the signal components to be included in 
the fits, the data were used as a guide. As the i<i(1270) 
signal is prominent in both B + — > J/iJjK + ir + ir~ and 
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TABLE V. Fitted parameters of the signal function for B + — > J/ipK + -iv + -K , along with the corresponding decay fractions. 



Jl 


Submode 




Modulus 


Phase (radians) 




Decay fraction 




Nonresonant K + n + n 


1.0 (fixed) 


(fixed) 


0.152 


± 


0.013 


± 


0.028 




#i(1270) - 


H> #*(892)tt 


0.962 ± 0.058 ± 0.176 


(fixed) 


0.232 


± 


0.017 


± 


0.058 




#i(1270) - 


#p 


1.813 ± 0.090 ± 0.243 


-0.764 ± 0.069 ± 0.127 


0.383 


± 


0.016 


± 


0.036 


1+ 


#i(1270) - 


#u; 


0.198 ± 0.036 ± 0.041 


1.09 ±0.18 ±0.18 


0.0045 


± 


0.0017 


± 


0.0014 




#i(1270) - 


-> #o(1430)tt 


0.95 ± 0.16 ± 0.24 


2.83 ±0.18 ±0.18 


0.0157 


± 


0.0052 


± 


0.0049 




#i(1400) - 


-> #*(892)tt 


0.894 ± 0.066 ± 0.125 


-2.300 ± 0.044 ± 0.078 


0.223 


± 


0.026 


± 


0.036 


1" 


A"* (1410) - 


^ #*(892)tt 


0.516 ± 0.090 ± 0.103 


(fixed) 


0.047 


± 


0.016 


± 


0.015 




# 2 *(1430) - 


-> #*(892)tt 


0.663 ± 0.051 ± 0.085 


(fixed) 


0.088 


± 


0.011 


± 


0.011 




# 2 *(1430) - 


^ #p 


0.371 (fixed) 


-1.12 ±0.22 ±0.29 




0.0233 (fixed) 




2+ 


# 2 *(1430) - 


-> #w 


0.040 (fixed) 


0.58 ±0.51 ±0.27 




0.00036 (fixed) 




if! (1980) - 


-> #*(892)tt 


0.775 ± 0.054 ± 0.118 


-1.59 ±0.15 ±0.14 


0.0739 


± 


0.0073 


± 


0.0095 




# 2 *(1980) - 


-> #p 


0.660 ± 0.048 ± 0.101 


0.86 ±0.22 ±0.21 


0.0613 


± 


0.0058 


± 


0.0059 




#(1600) - 


► #*(892)tt 


0.131 ± 0.021 ± 0.024 


(fixed) 


0.0187 


± 


0.0058 


± 


0.0050 




#(1600) - 


> #p 


0.193 ± 0.017 ± 0.029 


-0.27 ± 0.27 ± 0.18 


0.0424 


± 


0.0062 


± 


0.0110 


2" 


#2 (1770) - 


->■ #*(892)tt 


0.122 ± 0.021 ± 0.026 


2.22 ±0.49 ±0.37 


0.0164 


± 


0.0055 


± 


0.0061 




if 2 (1770) - 


->■ #|(1430)7r 


0.286 ± 0.043 ± 0.044 


1.78 ±0.39 ±0.24 


0.0100 


± 


0.0028 


± 


0.0020 




#2(1770) - 


-> #/ 2 (1270) 


0.444 ± 0.069 ± 0.077 


2.30 ±0.37 ±0.32 


0.0124 


± 


0.0033 


± 


0.0022 




#2(1770) - 


-7 #/o(980) 


0.113 ± 0.029 ± 0.024 


1.83 ±0.45 ±0.53 


0.0034 




0.0017 




0.0011 



TABLE VI. Fitted parameters of the signal function for B + — > tp 1 'K + ir + ir , along with the corresponding decay fractions. 



Jl 


Submode 




Modulus 


Phase (radians) 




Decay Fraction 




Nonresonant # + 7r + 7r 


1.0 (fixed) 


(fixed) 


0.253 


± 0.045 


± 


0.102 




#i(1270) - 


-> #*(892)tt 


0.213 ± 0.037 ± 0.049 


(fixed) 


0.090 


± 0.024 


± 


0.013 


1+ 


#1(1270) - 


■> Kp 


0.513 ± 0.070 ± 0.141 


-0.66 ± 0.26 ± 0.11 


0.215 


± 0.038 


± 


0.045 




#i(1270) - 


* #w 


0.048 ± 0.041 ± 0.022 


-0.37 ± 1.21 ± 0.52 


0.0017 


± 0.0033 


± 


0.0013 




#*(1680) - 


■» #*(892)tt 


0.67 ± 0.12 ± 0.15 


(fixed) 


0.106 


± 0.031 


± 


0.017 


1" 


#* (1680) - 


-> Kp 


1.17 ± 0.16 ± 0.32 


1.27 ± 0.24 ± 0.15 


0.241 


± 0.047 


± 


0.050 




#*(1680) - 


■» #CJ 


0.233 ± 0.097 ± 0.047 


-3.06 ± 0.43 ± 0.45 


0.0119 


± 0.0106 


± 


0.0061 



TABLE VII. Fitted parameters of the signal function for B + —¥ J/tpK + ir + -K when the #i(1270) mass and width are floated, 
along with the corresponding decay fractions. 



Jl 


Submode 




Modulus 


Phase (radians) 




Decay Fraction 






Nonresonant # + 7r + 7r 


1.0 (fixed) 


(fixed) 


0.142 


± 


0.013 


± 


0.026 




#i(1270) - 


-)■ #*(892)tt 


0.882 ± 0.076 ± 0.090 


(fixed) 


0.168 


± 


0.023 


± 


0.012 




#i(1270) -»■ Kp 


2.14 ±0.12 ±0.27 


-0.588 ± 0.084 ± 0.110 


0.430 


± 


0.018 


± 


0.027 


1+ 


#i(1270) - 


-> #CJ 


0.289 ± 0.043 ± 0.040 


1.25 ±0.16 ±0.14 


0.00758 


± 


0.00216 


± 


0.00076 




#i(1270) - 


-> #o(1430)tt 


1.09 ±0.18 ±0.24 


2.93 ±0.18 ±0.16 


0.0184 


± 


0.0055 


± 


0.0046 




#i(1400) - 


-> #*(892)tt 


0.746 ± 0.085 ± 0.089 


-2.585 ± 0.100 ± 0.076 


0.145 


± 


0.029 


± 


0.017 


\- 


#*(1410) - 


-> #*(892)tt 


0.736 ± 0.084 ± 0.098 


(fixed) 


0.089 


± 


0.019 


± 


0.010 




#2* (1430) - 


#*(892)tt 


0.529 ± 0.064 ± 0.070 


(fixed) 


0.0525 


± 


0.0120 


± 


0.0070 




# 2 * (1430) - 


-> Kp 


0.296 (fixed) 


-0.61 ±0.39 ± 1.07 




0.014 (fixe 


d) 




2+ 


#1(1430) - 


-> Ku) 


0.032 (fixed) 


1.41 ±0.80 ±0.25 




0.00021 (fixed) 






# 2 * (1980) - 


-> #*(892)tt 


0.756 ± 0.060 ± 0.119 


-1.46 ±0.23 ±0.22 


0.0659 


± 


0.0073 


± 


0.0088 




# 2 * (1980) - 


-> Kp 


0.685 ± 0.052 ± 0.106 


1.15 ±0.30 ±0.20 


0.0617 


± 


0.0061 


± 


0.0065 




#(1600) - 


> #*(892)tt 


0.147 ±0.021 ±0.026 


(fixed) 


0.0222 


± 


0.0058 


± 


0.0054 




#(1600) - 


> Kp 


0.171 ± 0.020 ± 0.023 


-0.21 ±0.29 ±0.12 


0.0312 


± 


0.0065 


± 


0.0040 


2" 


# 2 (1770) - 


#*(892)tt 


0.116 ± 0.020 ± 0.029 


1.93 ±0.52 ±0.43 


0.0137 


± 


0.0049 


± 


0.0058 




# 2 (1770) - 


-»■ # 2 *(1430)tt 


0.288 ± 0.045 ± 0.045 


1.81 ±0.40 ±0.24 


0.0095 


± 


0.0027 


± 


0.0018 




# 2 (1770) ->■ #/ 2 (1270) 


0.466 ± 0.073 ± 0.080 


2.22 ±0.36 ±0.33 


0.0128 


± 


0.0035 


± 


0.0021 




# 2 (1770) - 


->#/o(980) 


0.118 ± 0.030 ± 0.025 


1.89 ±0.45 ±0.53 


0.0035 


± 


0.0017 


± 


0.0011 
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0.60 GeV7c 4 < M 2 (K;t7i) < 1 .46 GeV7c 4 
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FIG. 22. B+ -> J/iPK+tt+tv- signal 
color coded as shown in Fig. 1191 
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0.60 GeV 2 /c 4 < M 2 (K;t7i) < 1 .00 GeV 2 /c 4 
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FIG. 23. B + — > ip'K + iv + TY signal data (points) and fit results (histograms) for slices in M 2 (Kttit). The fit components are 
color coded as shown in Fig. 1191 
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0.60 GeV 2 /c 4 < M 2 (K;m) < 1 .46 GeV 2 /c 4 
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FIG. 24. _B + — >■ J /ipK + tv + ty signal data (points) and fit results (histograms) for slices in M 2 (Kttw). The fit components are 
color coded as shown in Fig. I19I The mass and width of the _Ki(1270) floated in the fit. 
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B + —> t//K + iT + ir data, the initial fits were done with 
only #i(1270) -> A*(892)tt and ^(1270) -> A> on top 
of the nonresonant component. Additional decay chan- 
nels were added successively until a reasonable level of 
agreement between fit and data was obtained. 

As a further guide, the decays B° — > J/i/jK+tt^ 
and B° — > ip 1 ' K + tt~ were reconstructed. The observed 
Kit mass spectra are shown in Fig. [55] Consistent 
with the 1 + spin-parity assignment of the Ai(1270), no 
Ai(1270) — > Kn signal appears in these spectra. In 
both modes, a small peak can be seen near 1.4 GeV/c 2 
in M(Kw); this may have contributions from AT*(1410) 
or A 2 *(1430), as well as A *(1430). The absence of a 
A"* (1680) peak in B° -> J/^K + n~ is noteworthy, al- 
though a precise statement would require an analysis of 
the efficiency and phase space for these modeso In 
B° —> ip'K + Tr~, the kinematically-allowed M(Ktt) re- 
gion does not allow any conclusions to be drawn about 
the presence or absence of a low A* (1680) tail. 



Figure [27] shows scatterplots of signal-region B + —> 
J /ipK + ir + ix~ data over the three coordinates. Interfer- 
ence between Ai(1270) -> A*(892)tt and Ai(1270) -)• 
Kp is responsible for the weakening of the latter signal 
at M(Kir) > M K *( 8g2 )- Although the four-body phase 
space decreases with increasing M(Kir), this is not suf- 
ficient to account for the abrupt falloff. To describe the 
data in this region, the two modes must be added coher- 
ently. 

Since the previously-measured [Til ] branching fraction 
for Ai(1270) — > Kuj is small compared to that for 
Ai(1270) — > Kp, and since only 1.5% of uj's decay to 
7T + 7r |H3 one might expect Ai(1270) — > Kuj to play a 
negligible role in this analysis. Nonetheless, since the 
uj is much narrower than the p, it significantly distorts 
the observed p line shape through interference [27j ■ In 
Fig. OH the M 2 {-ktt) projections of Figs. OH [20] and ED 
are finely binned to demonstrate this interference pat- 
tern, which is accurately modeled by the PDFs. 
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FIG. 25. Observed Kit mass spectra for B° 
(top) and B° -> tp'K + -R~ (bottom) data. 
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2. Interference effects 

The inclusion of interference among submodes sharing 
the same initial-state spin-parity is essential to obtaining 
good fits to the data. In particular, dramatic interference 
effects are observed between Ai(1270) -)• A*(892)7r and 
Ai(1270) -> Kp, as well as between Ai(1270) -> Kp and 
Ai(1270) -)■ Kuj. 
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A detailed analysis of B — ¥ J /ipK-K and B — > ip' Kit is beyond 
the scope of this work. A Dalitz analysis of the latter mode was 
presented in Ref. l25i 



Although u decays dominantly to to it + tt~ it , it can also decay 
to tt + tt~ through G-parity violation |26| , which causes mixing 
between p and uj. An u> component is therefore present whenever 
a particle decays to p. 
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FIG. 26. Finely-binned projections onto the M 2 (tttt) axis 
of signal-region data (fits) and fit results (histograms) for 
B+ -» J/tpK+TT+TT- (top), B+ -> i/''-K' + 7r + 7r" (middle), and 
B+ -> J/TpK + TT + TT~ with the mass and width of the Id (1270) 
floated (bottom). The fit components are color coded as 
shown in Fig. 1131 The discontinuity at the u mass is due 
to p-u) interference. 

The peculiar shape of the observed p-uj interference 
pattern is caused by kinematic effects. The largest con- 
tribution to the p signal comes from K i (1270) — > Kp, 
which straddles the edge of phase space, as can be seen in 
the middle panel of Fig. [37] The distortion that is caused 
by this kinematic cutoff is taken into account automat- 
ically by integrating the signal function only over the 
kinematically-allowed region, as described in Sec. IVI Al 
Modeling the data accurately requires including p-uj in- 
terference, incorporating the four-body phase space fac- 
tor into the signal function, and integrating the signal 
function over only the kinematically-allowed phase space. 



3. The L region 

The M(Kmr) region between 1.5 and 2.0 GeV/c 2 , his- 
torically referred to as the L reg ion, comprises several 
wide, overlapping resonances [28j~30]. The large uncer- 
tainties in the masses and widths of the known states in 
this region make it difficult to characterize this region in 
this analysis. The model presented here is not necessarily 
the only one supported by the data. 

To describe the structure observed at 2.6 GeV 2 /c 4 
in the M 2 (Ktttt) distribution of B + -)• J/^jK+tt+tt^ , 
a peak with a mass of 1.605 GeV/c 2 and a width of 
115 MeV/c 2 is included in the fit, decaying to #*(892)tt 
and Kp. This peak, which is referred to as if (1600) in 
this paper, may be the #2(1580), an as-yet unconfirmed 
J p = 2~ state that has previously been observed decay- 
ing to #* (892)tt [H[. 

As can be seen in Fig. [22] the high end of the M 2 (Ktttt) 
spectrum of B + — > J/i(jK + tt + tt~ exhibits #* and p sig- 
nals. To fit the data in this region, we include a #2 (1980) 
resonance, which is another state that currently requires 
confirmation. 

Even after including #(1600) and if* (1980) reso- 
nances in the B + — > J/ipK + ir + ir~ fit, a slight en- 
hancement remains around 3 GeV 2 /c 4 in M 2 (K TTTT). 
A #2(1770) signal is therefore also included, with its 
known decays to #*(892)tt, # 2 *(1430)tt, #/ (980), and 
#/ 2 (1270). 

Fitting the B + — > tp' K + tt + tt~ data is more difficult 
still, as there are fewer events to analyze, and only a 
small portion of the L-region is within the kinematic 
limits of the decay. In addition to the #i(1270) signal, 
the M 2 (Ktttt) spectrum contains what appears to be the 
low-mass tail of at least one high-mass resonance. As 
Fig. [53] shows, there are clear K* (892) and p peaks at 
high M 2 (Ktttt); these are not reproduced by the PDF 



if no high-mass resonance is included in the model. If 
the enhancement is modeled as a single resonance, the 
data favor a mass of roughly 1.7 GeV/c 2 and a width 
of 400-500 MeV/c 2 . In this analysis, the enhancement 
is modeled as the #*(1680). The data do not preclude 
other possibilities, such as the #2(1770). Indeed, the hint 
of / (980) in the last M 2 (#7T7t) slice in Fig. |2"31 cannot 
come from a 1~ state such as the #*(1680), or from a 
2+ state such as the #2 (1430). 



4- Comparison with previous measurements 

It is interesting to compare the relative decay fractions 
for #i(1270) submodes in the B+ J/^K+tt+tt- fits to 
previous measurements of #i(1270) branching fractions. 
For this purpose, we use the decay fractions with phase 
space shown in Tables [V] and I VIII include isospin fac- 
tors, and assume branching fractions of (l-53lg_j_i)% for 
lo ->• 7T+7T-, and (93±10)% for # *(1430) -» Ktt Q. The 
calculation neglects the systematic errors in the decay 
fractions and assumes that the statistical errors among 
the decay fractions are uncorrelated. Moreover, it as- 
sumes that the #i(1270) decays only to #*(892)tt, Kp, 
Koj, and #o(1430)-7r, and neglects interference among 
these decay channels. The comparison is shown in Ta- 
ble ELLD While the ratios of the #i(1270) branching 
fractions to #*(892)7r, Kp, and Kui are consistent with 
the previously-measured values, the branching fraction 
to #g(1430)7T is significantly smaller. 



TABLE VIII. Comparison of branching fractions for #i(1270) 
decays according to the Particle Data Group (PDG) [lj] and 
based on the results shown in Table IVl (fit 1) and Table fVTTl 
(fit 2). See text for assumptions. 



#i(1270) 








Branching Fraction (%) 


Decay mode 


PDG 


Fit 1 


Fit 2 


Kp 


42 


± 


6 


57.3 ± 3.5 


58.4 ± 4.3 


# *(1430)7T 


28 


± 


4 


1.90 ± 0.66 


2.01 ± 0.64 


#*(892)tt 


16 


± 


5 


26.0 ± 2.1 


17.1 ± 2.3 


Koj 


11 


± 


2 


14.8 ± 4.7 


22.5 ± 5.2 


#/o(1370) 


3 


± 


2 


N/A 


N/A 



5. Mass and width of the A"i(1270) 

As shown in Sec. IVI II the data favor a smaller mass 
and a larger width for the #i(1270) than the Particle 
Data Group (PDG) values. This is mainly due to the 
excess of #*(892) and p at low M 2 (Ktttt), as can be 
ascertained by comparing the first row of plots in Figs. [22] 
and 1241 The measured mass and width agree remarkably 
well with Ref. [3l| and are also consistent with Ref. HH. 
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FIG. 27. Scatterplots for B H 



J/ipK + Tv + -K data, showing M 2 (iT7r) versus M 2 (Kmr) (left), M 2 (7T7r) versus M 2 (K 7T7T 



(middle), and M 2 (Kir) versus M 2 (nir) (right). Interference between the Kp and if*(892)-7r submodes of the ifi(1270) is 
responsible for the abrupt fading of the if* (892) signal at M(K-k) > M K *( 8 g 2 )- The effect is most apparent in the left and 
right plots. Scatterplots for B + —¥ ip'K + -K + ir~ data are similar but limited by statistics. 



6. Limitations of the method 

There are large uncertainties in the masses and widths 
of many of the states included in the fits, as can be seen 
in Table IIV1 Although this is taken into account in cal- 
culating the systematic error, it nonetheless limits the 
accuracy of the model. 

In B + — > ip' K + ir + Tr~ , the small sample size and the 
kinematic cutoff limit the conclusions that can be drawn 
about the signal components. In B + — > J/ipK + ir + ir~ , 
the sample size is larger, but a further limitation is im- 
posed by the increase in computation time as more pa- 
rameters are added to the fit. Each additional decay 
channel that is included in the signal function contributes 
a modulus and possibly a phase to be varied in the fit. 
Since the normalization integral of the signal function in 
Eq. [TC] depends on the values of the parameters a, the 
integration must be performed for each set of parame- 
ters attempted by the fitter. While the step size used in 
the numerical integration can be increased to speed up 
the process, it must be small enough to allow the PDF 
to resolve the structures in the data. In particular, the 
p-uj interference pattern can be fitted with a step size of 
0.01 GeV 2 /c 4 , but not with a step size of 0.02 GeV 2 /c 4 . 
As a consequence of the finite processor speed, not every 
possible decay channel can be included in the fit. The 
model is necessarily incomplete. 

The large nonresonant component seen in both B + 
J/ipK + 'K + Ti~ and B + — >• ip 1 K + Tr + Tt~ may be an indi- 
cation of contributions from additional wide kaon exci- 
tations. It may also incorporate some misreconstructed 
resonant signal. While the nonresonant component is 
assumed in this analysis to be distributed according to 
phase space, this assumption may be inaccurate. There 
are currently no accepted models of nonresonant B- 
meson decays. 

It is difficult, in an analysis like the one presented here, 
to determine the significance of a given component of the 



signal. An improvement in the likelihood upon the addi- 
tion of a new resonance to the signal function indicates 
only that the model is incomplete, not necessarily that 
the data contain the particular resonance. Furthermore, 
unless the model is accurate in every other way, floating 
the mass and width of a particle in the fit may not yield 
a reliable result, as the fitter may set these parameters 
to compensate for the model's deficiencies. This is espe- 
cially important in the high-M 2 (A'7r7r) region, where the 
statistics are limited and there are large uncertainties in 
the masses and widths of the resonances included in the 
signal function. Thus, although the if (1600) component 
of the signal function for B + — > J/ipK + Tr + ir~ greatly im- 
proves the quality of the fit, it is difficult to claim that it 
is a single particle, let alone measure its mass and width. 



VII. CONCLUSIONS 

Using data recorded by the Belle detector, we have 
measured branching fractions for the decays B + — > 
J /ipK + K + n~ and B + — > ip' K + Tr + ir^ with improved pre- 
cision (see Sec. IV Bp . We have also performed ampli- 
tude analyses in three dimensions — M 2 (Kirn) , M 2 (Kir) , 
and M 2 (tttt) — to determine the resonant structure of the 
K + tt + tt~ final state in these decays (see Sec. IVI II) . 

We have shown that the ifi(1270), which is the dom- 
inant component of the K + -k + tt~ final state in B + — > 
J/tpK + ir + ir~ , is also prominent in B + — > ip' K + ir + ir~ . 
The large sample available for the former decay reveals 
a small peak at M(Ktttt) « 1.4 GeV/c 2 . Our three- 
dimensional fits represent a first attempt to determine 
the components of this peak. 

Performing an unbinned fit in three dimensions ex- 
ploits practically all of the information available in the 
data. While it is relatively easy to obtain a good fit in 
one dimension, requiring a fit that succeeds in three di- 
mensions greatly restricts the class of successful models. 
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With high statistics, it is possible to use interference ef- 
fects and the spin-dependent angular distribution of the 
final state to distinguish overlapping resonances. In par- 
ticular, we have shown that p-uj interference cannot be 
neglected in studying (1270) decays to Kir + ir~ final 
states. 

The large size of the B + — > J /ij)K + iT + ir~ data sample 
allows us to measure the mass and width of the i£"i(1270) 
with improved precision (see Eqs. [38] and [39]). These 
values differ considerably from previously-published val- 
ues |14|. The analysis of these data also provides infor- 



mation on the relative strengths of A"i(1270) decays to 
Kp Kw , iT*(892)7r, and JC *(1430)7r final states (see Ta- 
bic [VIlI|. While the results are consistent with previous 
measurements for the first three modes, they indicate a 
much smaller rate of decay to Kq (1430)tt than previously 
accepted. 

Although more data are required to clarify the struc- 
ture of the high M 2 (K + tt + it~) region in both B + — > 
J/ipK + Tt + Tr~ and B + — > ip' K + tt + t:~ , we have shown 
that this region contains broad resonances that decay to 
A*(892)tt and Kp final states. 

The analysis presented in this paper demonstrates 
that the decay modes B + — > J/ijjK + ir + ir~ and B + — > 
■ij/ K + ir + Tr~ can provide clean laboratories for the spec- 
troscopy of excited kaon states. Many of these states 
still require confirmation or more precise mass and width 
measurements. As more data become available at future 
super-B factories, analyses similar to the one presented 
here can further elucidate the higher regions of the kaon 
spectrum. 
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APPENDIX: LOCAL MAXIMA 

This appendix summarizes the results of the local- 
maximum test, in which each of the three signal-region 
fits was repeated 100 times with randomly selected start- 
ing values for the parameters. In each case, the best like- 
lihood obtained coincided with the solution presented in 
Tables [VTlVIH In the following, these solutions are re- 
ferred to as the "global maxima." 

For the B+ -> J/i/jK+tt+tt- fit with the # x (1270) 
mass and width fixed to their values in Table IIV1 the 
local maximum closest to the global maximum presented 
in Table [V] had a likelihood of -10608.6, which is 8.2cr 
away from the global maximum. 

For the B + — s> ip' K + tt + tt~ fit, two local maxima were 
found: one with a likelihood of 635.4 and the other with 
a likelihood of 635.9; these are 2.4<r and 2.2<r away from 
the global maximum, respectively. The former had an 
unphysically large decay fraction for iTi(1270) — > Klj 
and was discarded. For the latter, all the parameters 
were within statistical error of the values presented in Ta- 
ble EH with the exception of the Ai(1270) -> A*(892)tt 
amplitude and decay fraction, which were higher by 1.4 
times the statistical error. 

For the B+ -» J/iPK+n+ir- fit with the Ai(1270) 
mass and width allowed to float, two local maxima were 
found, both with a likelihood of —10528.8, which is 2.7<7 
away from the global maximum. The fitted parameters 
and decay fractions for these local maxima are presented 
in Tables IIXI and [X[ respectively. The fitted mass and 
width are 

M Kl(1270) = (1241.9 ±3.2) MeV/c 2 , 
r^cmo) = (128.3 ±5.8) MeV/c 2 



for the former, and 



M 



K"i(1270) 



K i (1270) 



(1244.3 ±3.3) MeV/c 2 , 
= (129.0 ±5.7) MeV/c 2 



for the latter. 
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TABLE IX. Fitted parameters and decay fractions corresponding to the first local maximum for the mode B + — ¥ J/ipK + -K + ir , 
with the if 1 (1270) mass and width floated. The errors are statistical. 



Jl 


Submode 




Modulus 


Phase (radians) 


Decay Fraction 




Nonresonant K + n + n 


1.0 (fixed) 


(fixed) 


0.139 ± 0.015 




if i (1270) -S- if*(892)7r 


0.662 ± 0.056 


(fixed) 


0.090 ± 0.012 




ifi(1270) - 


-> if p 


2.22 ± 0.13 


-0.96 ± 0.13 


0.438 ± 0.021 


1+ 


if i (1270) - 


-> ifw 


0.301 ± 0.043 


0.94 ± 0.18 


0.0078 ± 0.0022 




if i (1270) - 


^(1430)71 


0.88 ± 0.18 


2.39 ± 0.27 


0.0116 ± 0.0047 




if i (1400) - 


-> if*(892)7r 


0.258 ± 0.083 


2.35 ± 0.43 


0.017 ± 0.011 


1" 


A"* (1410) - 


-> if*(892)7r 


0.755 ± 0.099 


(fixed) 


0.091 ± 0.023 




if 2 * (1430) - 


if*(892)7r 


0.384 ± 0.081 


(fixed) 


0.027 ± 0.012 




if 2 *(1430) - 


+ ifp 


0.214 (fixed) 


2.63 ± 0.60 


0.0071 (fixed) 


2+ 


if 2 * (1430) - 


-> ifw 


0.023 (fixed) 


-2.6± 1.1 


0.00011 (fixed) 




if| (1980) - 


-> if*(892)7r 


0.659 ± 0.059 


0.42 ± 0.48 


0.0487 ± 0.0064 




^5(1980) - 


-> ifp 


0.733 ± 0.057 


2.92 ± 0.40 


0.0689 ± 0.0071 




if (1600) - 


► if*(892)7r 


0.175 ± 0.066 


(fixed) 


0.031 ± 0.021 




if (1600) - 


► ifp 


0.169 ± 0.020 


0.61 ± 0.50 


0.0297 ± 0.0070 


2" 


if 2 (1770) - 


-> if*(892)7r 


0.138 ± 0.066 


2.9± 1.3 


0.019 ± 0.017 




if 2 (1770) -S- if 2 *(1430)7r 


0.305 ± 0.055 


2.55 ± 0.68 


0.0104 ± 0.0032 




if 2 (1770) - 


-S- Kj 2 (1270) 


0.515 ± 0.090 


2.81 ± 0.79 


0.0153 ± 0.0045 




if 2 (1770) - 


-*if/o(980) 


0.121 ± 0.031 


2.54 ± 0.69 


0.0036 ± 0.0019 



TABLE X. Fitted parameters and decay fractions corresponding to the second local maximum for the mode B + — > J/ipK + -K + ir , 
with the if i (1270) mass and width floated. The errors are statistical. 



Jl 


Submode 




Modulus 


Phase (radians) 


Decay Fraction 




Nonresonant K + n + n 


1.0 (fixed) 


(fixed) 


0.136 ± 0.013 




ifi(1270) - 


-S- if*(892)7r 


0.701 ± 0.060 


(fixed) 


0.100 ± 0.013 




if i (1270) ifp 


2.24 ± 0.13 


-1.00 ± 0.10 


0.447 ± 0.018 


1+ 


if i (1270) - 


-> Kuj 


0.293 ± 0.043 


0.91 ± 0.16 


0.0075 ± 0.0020 




if i (1270) - 


-> ^(1430)71 


0.91 ± 0.18 


2.40 ± 0.25 


0.0123 ± 0.0045 




if i (1400) - 


if*(892)7r 


0.162 ± 0.076 


2.53 ± 0.45 


0.0066 ± 0.0060 


1" 


if* (1410) - 


-> if*(892)7r 


0.739 ± 0.093 


(fixed) 


0.086 ± 0.020 




if 2 * (1430) - 


-+ if*(892)7r 


0.687 ± 0.082 


(fixed) 


0.085 ± 0.018 




if 2 *(1430) - 


-> Kp 


0.384 (fixed) 


3.01 ± 0.27 


0.022 (fixed) 


2+ 


if 2 * (1430) - 


Kuj 


0.041 (fixed) 


-1.4± 1.2 


0.00034 (fixed) 




if 2 * (1980) - 


-> K*(892)n 


0.814 ± 0.063 


-1.52 ± 0.19 


0.0731 ± 0.0077 




iff (1980) - 


Kp 


0.791 ± 0.059 


0.81 ± 0.25 


0.0788 ± 0.0070 




if (1600) - 


► if*(892)7r 


0.253 ± 0.043 


(fixed) 


0.063 ± 0.019 




if (1600) - 


► Kp 


0.152 ± 0.021 


0.86 ± 0.29 


0.0236 ± 0.0062 


2" 


if 2 (1770) - 


if*(892)7r 


0.221 ± 0.050 


-2.56 ± 0.26 


0.048 ± 0.020 




if 2 (1770) - 


-7 if 2 *(1430)7r 


0.326 ± 0.048 


2.98 ± 0.40 


0.0117 ± 0.0030 




if 2 (1770) if/ 2 (1270) 


0.541 ± 0.073 


-2.89 ± 0.34 


0.0166 ± 0.0037 




if 2 (1770) -> if/o(980) 


0.121 ± 0.032 


2.98 ± 0.48 


0.0035 ± 0.0018 
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